Tulai

Author

S.SH

last updated

July 18, 2023

# library
library(tidyverse)
library(ggplot2)
library(dplyr)
library(ggbeeswarm)
# raed the data

tl <- read_csv("Tulai_Lithic_Assamblage.csv")
# create a new data frame for plotting
tl1 <- tl %>%  
  select(AREA,
         `Width (mm)`, 
         `Thickness (mm)`,
         `Length (mm)`)  
tl1 %>%
  filter(AREA == "A") %>%
  select(AREA, `Width (mm)`) %>%
  ggplot() +
  aes(x = `Width (mm)`) +
  geom_histogram(binwidth = 2) +
  labs(x = "Width of Debitage (mm)", 
       y = "Frequency", 
       title = "Histogram of Debitage Width for Area A" )+
  theme_minimal()

Figure 1: Width of Artifacts (Area A)
#| label: Width of Artifacts (Area A)
#| fig-cap: "Width of Artifacts (Area A)"
#| fig-height: 5

tl1 %>%
  filter(AREA == "TP1") %>%
  select(AREA, `Width (mm)`) %>%
  ggplot() +
  aes(x = `Width (mm)`) +
  geom_histogram(binwidth = 2) +
  labs(x = "Width of Debitage (mm)", 
       y = "Frequency", 
       title = "Histogram of Width for Area TP1" )+
  theme_minimal()

tl1 %>%
 filter(AREA == "D1, D3, D4, D5, C, B1, DH") %>%
  select(AREA, `Width (mm)`) %>%
  ggplot() +
  aes(x = `Width (mm)`) +
  geom_histogram(binwidth = 5) +
  labs(x = "Width of Debitage (mm)", 
       y = "Frequency", 
       title = "Histogram of Width for other Areas" )+
  theme_minimal()

ggplot(data = tl1) + 
       aes(x = `Width (mm)`,  
           y = AREA) +  
  geom_boxplot() +  
  labs(x = "Width",  
       y =  "Areas") +  
  theme_classic(base_size = 14)